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I  • 


1 .  INTRODUCTION 


Since  the  publication  of  Anderson's  classic  paper, *  questions 
concerning  impact  theories  of  pressure  broadening  in  the  microwave  and 
infrared  regions  have  shifted  from  the  realm  of  the  foundations  of  the 
impact  theory.  Attention  is  now  concentrated  on  such  areas  as  the 
convergence  of  the  various  terms  in  Anderson’s  expansion  and  the 
extension  of  the  theory  to  include  broadening  of  vibration-rotation 
lines,  as  well  as  the  pure  rotation  lines  addressed  by  Anderson.  Ibe 
solutions  of  these  problems,  particularly  the  former,  have  not  been 
entirely  satisfactory.  Several  ad  hoc  cutoffs  have  been  introduced  into 
Andersons’s  theory, 1-3  and  although  these  have  agreed  well  with  each 
other  and  with  experiment,  their  differences  are  sufficient  to  warrant 
further  examination. 

In  an  attempt  to  resolve  these  problems,  Murphy  and  Boggs4  have 
developed  a  theory  which  begins  with  a  different  theoretical  framework 
from  Anderson's,  and  which  introduces  a  natural  cutoff  from  the 
outset.  This  theory  has  several  drawbacks,  however.  First,  the  cutoff 
depends  on  the  azimuthal  quantum  numbers  of  the  radiating  and  perturbing 
molecules.  This  requires  an  involved  machine  computation  to  evaluate 
certain  m-sums  occurring  in  the  theory.  Second,  the  theoretical 

foundation  of  the  theory  is  not  as  clear  as  that  of  Anderson;  no  off- 
diagonal  matrix  elements  of  the  time-evolution  operator  occur  in  Murphy 
and  Boggs'  theory,  and  it  is  not  clear  how  these  are  to  be  introduced 
into  their  framework. 

*P.  W.  Anderson,  Phys,  Rev.,  76_  (1949),  647. 

2P.  W.  Anderson,  PhD  Dissertation,  Harvard  University  (1949). 

3P.  D.  Sharma  and  G.  E.  Caledonia,  J.  Chem.  Phys.,  _5 £  (1970) ,  434. 

4J.  S.  Murphy  and  J.  E.  Boggs,  J.  Chem.  Phys.,  47  (1967),  691. 
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Gordon  has  given  a  nonperturbative  semiclassical  theory#5  which 
includes  all  the  important  effects  of  molecular  collisions:  phase 
shifts,  inelastic  collisions,  and  molecular  reorientation  (the  Murphy- 
Boggs  theory  does  not  include  the  third  of  these  effects).  However,  in 
order  to  evaluate  pressure-broadening  cross  sections  from  his  theory, 
Gordon  is  forced  into  a  rough  guess  for  the  probability  of  elastic 
collisions  (because  he  treats  the  angular  momentum  as  a  continuous 
variable).  Thus,  Gordon’s  theory,  while  satisfactory  from  a  theoretical 
point  of  view,  does  not  allow  accurate  calculations  of  pressure  shifts 
and  broadening. 

What  we  propose  is  an  impact  theory  of  pressure  broadening  and 
shift,  based  on  Anderson's  framework,  and  containing  a  natural,  m- 
independent  cutoff.  The  theory  is  similar  to  that  of  Murphy  and  Boggs, 
except  that  it  explicitly  exhibits  the  correct  rotational  symmetry.  The 
theory  is  expressed  from  the  outset  in  terms  of  the  group-theoretical 
properties  of  the  interaction  Hamiltonian  between  the  radiating  and 
perturbing  molecules.  In  this  way  it  is  possible  to  formulate  the 
theory  in  a  very  general  manner,  allowing  for  complicated  interactions 
between  the  two  molecules.  The  final  form  of  the  theory  is  shown  to  be 
similar  to  the  form  of  the  semiclassical  theory  of  Gordon. 

Our  starting  point  will  be  the  general  relations  given  by  Anderson1 
and  rederived  by  Tsao  and  Curnutte.5  In  particular,  we  use  equations 
(79)  and  (80)  of  Tsao  and  Curnutte : 

nv  -l 

Av  *  -  o,  cm  1  ( 1 ) 

2ttc  1 


lP.  ff.  Anderson,  Phys,  Rev.,  76_  (1949),  647. 

5R.  G.  Gordon,  J.  Chem.  Phys.,  44_  (1966),  3083. 

6C.  J.  Tsao  and  B.  Curnutte,  J.  Quant.  Spect.  Rad.  Trans.,  2  (1962), 


41 


for  the  line  shift,  and 


(Av) 


1/2 


-  c_  cm 

2  ire  r 


(2) 


for  the  pressure-broadened  half-width.  In  the  above,  n  is  the  perturber 
density  and  v  is  the  relative  velocity  between  colliding  molecules.  The 
question  of  averaging  equations  (1)  and  (2)  over  a  Maxwell-Boltzmann 
distribution  of  relative  velocities  will  be  addressed  later.  For  now, 
we  simply  take  v  as  the  average  velocity. 


Anderson  gives  the  following  expression  for  the  cross  section: 


where 


a 


a  + 
r 


P 


J2  J2 


“2 


2ir 


b  db  S,  (b) 
J2 


(3) 


(4) 


and  where  J 2  is  the  rotational  angular  momentum  of  the  perturbing 
molecule,  and  b  is  the  impact  parameter.  The  latter  two  equations  are 
equations  (74)  and  (87)  of  Tsao  and  Curnutte.6  The  quantity  Pj2j2  a 
diagonal  element  of  the  density  matrix  for  the  perturber,  given  by 


6C.  J. 


41. 


Tsao  and  B. 


Curnutte ,  J.  Quant. 


Spect.  Rad.  Trans.,  2_  (1962) 
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where  ej2  is  the  energy  of  the  perturber  rotational  state  J2«  The 
function  Sj2(b)  is  given  by  equation  (88)  of  Tsao  and  Curnutte: 


S(b)  =  1  - 


v 

)1(M) u 

'l(Mi)>< 

(KM)  1  Ji(M[)> 

all  M  1 

[2Ja  +  1)(2J2  +  1] 

1 

(6) 


x  <JfMf  J2M2|T-1(b)  <JtM£J^|T(b)  |JiMiJ2M2> 


In  the  above,  the  subscripts  i  and  f  refer  to  initial  and  final  states 
of  the  radiative  transition,  and  primed  and  unprimed  quantities  refer  to 
states  after  and  before  the  collision,  respectively.  The  quantity 
^Jf  (Mf)  1  (M)  I  Ji  (Mi)  y  ■‘•s  an  unsymmetrized  vector-coupling  coefficient 
( Clebsch-Gordan  coef f icient) ,  and  T(b)  is  the  time-evolution  operator 
for  collisions  with  an  impact  parameter  b.  This  operator  obeys  the 
differential  equation 


ifi 


dT  ( b ;  t ) 
dt 


V(b;t)  T( b; t) 


subject  to  the  initial  condition 


(7) 


T( b; -“)  =  1 


(8) 


and 


T(b)  =  T(b,«) 


(9) 
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V(t)  is  the  interaction  Hamiltonian  between  radiator  and  perturber  in 
the  interaction  representation 


V(b;t)  =  exp  [-Hgt/i-rt]  V(t)  exp  [HQt/i-fi]  ,  (10) 


where  HQ  is  the  unperturbed  Hamiltonian  Jhq ( 1 )  +  HQ ( 2  fj  of  radiator  and 
perturber.  The  operator  V(t)  is  given  in  terms  of  operators  depending 
on  molecules  1  and  2  (radiator  and  perturber,  respectively),  and 
contains  an  explicit  time  dependence  due  to  the  changing  distance  and 
relative  orientation  between  the  two  molecules.  This  is  the  classical 
path  assumption.  We  have  suppressed  the  b-dependence  in  V(t);  such 
dependence  will  be  understood  throughout  the  following. 

We  take  equations  (1)  through  (10)  as  our  starting  points,  as  does 
Anderson's  theory.  However,  Anderson  goes  one  step  further  and  solves 
(7)  via  the  usual  perturbation  method.  He  obtains 

<m|T(b)  |n>  =  6^  +  ^—  I  dti  (mlV^)  |n> 

J-09 

(ID 

2  r°°  ft  i 

+  (i«)  /  dtl  /  at2  l  <»|v(t1)|k><t|v(t2)  ln>  +  .  .  • 

'  J—OO  J  —00  k 

where  we  have  used  the  shorthand  notation  (m)  to  stand  for  the  states 

(J1M1J2M2>.  Anderson  substitutes  expression  (11)  into  equation  (6)  to 

obtain  his  final  result  for  S„  (b).  When  this  expression  is  used,  the 

J2 

integral  of  equation  (4)  over  impact  parameter  diverges  due  to  the  small 
b  behavior  of  the  T-matrix  elements. 
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In  order  to  remove  this  divergence  in  a  natural  manner  (i.e., 
without  introducing  any  additional  ad  hoc  assumptions),  we  shall  obtain 
an  alternative  expansion  to  (11)  for  the  matrix  elements  of  the  T(b) 
operator.  To  do  this  we  examine  the  rotational  and  time-reversal 
symmetries  of  the  interaction  Hamiltonian  and  the  consequences  of  these 
symmetries  for  a  T-matrix  element  with  an  exponential  cutoff  factored 
out.  By  matching  this  exponential  expansion  with  the  expansion  of 
equation  (11),  we  can  unambiguously  determine  all  the  terms  in  the 
expansion.  We  also  show  that  the  exponential  expansion  is  equivalent  to 
summing  all  closed-loop  graphs  in  a  graphical  representation  of  the 
theory  by  use  of  a  linked  cluster  theorem. 


2.  GENERAL  DERIVATION 


2.1  Properties  of  T-matrix 

We  begin  our  derivation  of  the  exponential  cutoff  theory  by 
postulating  the  following  form  for  the  T-matrix  elements. 


<J|MiJ2M2|T|JtMlJ2H2>  = 


=  A 


,J1M1J2M2;J1M  1J2M2 


(12) 


1J2 


M2] 


where  the  quantities  A  and  B  are  yet  to  be  determined.  As  in  the 
introduction,  we  explicitly  label  the  states  according  to  their 
rotational  properties;  vibrational  and  electronic  quantum  numbers  will 
be  suppressed  in  what  follows.  We  also  assume  that  A  and  B  are  expanded 
in  a  power  series  in  the  interaction  potential.  If  the  right-hand  side 


t 
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of  equation  (12)  were  expanded  out  and  matched  term  by  term  with 

expression  (11),  a  system  of  equations  would  result  for  the  A^n'and 

jj(n)  [A^n^  represents  the  nt'^1  tern  in  the  power-series  expansion  for 

A,  and  similarly  for  B^nM«  However,  the  first  n  of  these  equations 

( k)  ( k) 

would  contain  2n  unknowns,-  that  is,  A  and  B  for  k  =  0,1,  ...,n-  . 

1.  Therefore,  some  other  conditions  on  the  A^n*  and  B^n^  are  necessary 
to  determine  them  uniquely. 

First  we  note  that  equation  (12)  must  be  substituted  into 
equation  (6)  to  obtain  S(b),  and  summations  must  be  carried  out  over  the 
M  quantum  numbers.  To  do  this  in  closed  form,  the  M  dependence  in  (12) 
must  be  of  the  form  of  vector-coupling  coefficients  or  products  of 
vector-coupling  coefficients,  and  these  must  occur  in  A  and  not  in  B. 
If  this  is  the  case,  we  may  use  the  well-known  properties  of  the  V-C 
coefficients  to  carry  out  the  M-summations .  Therefore,  we  postulate  the 
relation 


“  BJ' J ' ;  J  J  only 
1  1  2  2’  1  1  2  2  1212 


(13) 


This  relation  implies  that  the  quantity  j  is  rotationally 

invariant;  in  the  language  of  group  theory,  it  is  constructed  from 
combinations  of  the  interaction  Hamiltonian  which  transform  like 
irreducible  tensors  of  rank  zero  under  rotations  of  molecule  1  or 
molecule  2. 


For  our  second  relation,  we  examine  the  time-reversal  symmetry 
of  the  interaction  Hamiltonian.  We  assume  that  the  interaction 
Hamiltonian  may  be  expanded  in  terms  of  spherical  harmonics  in  the 
sector  R(t)  between  the  two  molecules. 


(14) 


* 


V(t)  =  l  AL/Pt(i,2)  yLy[e(t) ,4>(t)]  fL[R(t)] 
I*»u 


where  AL,ll(1,2)  is  a  time- independent  operator. 

To  study  the  time-reversal  properties  of  (14),  we  must  make 
some  assumptions  about  the  path  of  molecule  2.  It  is  usually  assumed 
that  the  molecule  travels  in  a  straight  line  relative  to  molecule  1; 
this  is  usually  a  good  approximation.  However,  we  need  not  be  so 
stringent,  since  we  wish  to  keep  the  theory  as  general  as  possible.  We 
shall  assume,  however,  that  the  translational  motion  of  the  two 
molecules  is  governed  by  central  forces,  and  that  the  noncentral  part  of 
the  potential  will  have  negligible  effect  on  the  translational  motion. 
If  this  is  the  case,  then  the  orbit  of  molecule  2  relative  to  molecule  1 
will  be  symmetric  about  the  point  of  closest  approach,  which  we  label  by 
the  vector  r0.  To  demonstrate  the  effect  of  this  orbit  symmetry  on  the 
symmetry  of  the  interaction  Hamiltonian  (14),  we  construct  a  specific 
coordinate  system  in  which  the  z-axis  is  coincident  with  the 
vector  rg  and  the  motion  of  molecule  2  is  in  the  y-z  plane.  We  also 
take  the  vector  R( 0 )  =  rg  so  that  the  orbit  is  symmetric  about  the  time 
t  =  0  in  the  following  way. 


R(-t)  =  R(t) 


(  15a) 


0(-t)  =  9(t) 


(15b) 


$(-t)  =  $(t)  +  ir 


(15c) 
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v  i 


where  0(t)  and  $(t)  are  the  usual  polar  angles  of  the  vector  R(t)  with 
respect  to  the  coordinate  system.  For  this  special  choice  of  coordinate 
system,  we  have 


=  YLU[0(t)'  *(t>] 


(16) 


Therefore 


hv2”2>* 

(17) 

=  <JiMiJ2M2lV(t)|JlMlJ2M2>  * 

To  see  the  consequences  of  this  relation,  we  examine  the  perturbation 
expansion,  (11).  By  reversing  the  order  of  integrations  and  by  using 
the  hermiticity  relation, 

<n|V(t) |m>  =  ln>*  >  (18) 

and  by  changing  the  signs  of  all  the  integration  variables,  we  arrive  at 
the  relation 

<n|T|m>=  6^  +  ^  /  dtj  <m  |V  (“tj)  |n>*  + 

a> 

(19) 

+  (in)2  /  dtl  /  1  dt2l  <"«|V(-ti)  |k>*<k|V(-t2)  |n>*  +  ... 

oo  •/■a)  lr 
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Substitution  of  (17)  into  (19)  yields  the  relation 


<JlM[J2M2lT|J1M1J2M2>  =  <J1M1J2M2!TlJiMiJ2M2> 


(20) 


Referring  to  equation  (12),  we  demand  that  (20)  hold  separately  for  each 
of  the  quantities  A  and  B.  The  relation  (20)  holds  for  T  only  in  the 
specific  coordinate  system  described  above.  However,  since  B  is 
rotationally  invariant,  the  symmetry  (20)  will  hold  for  B  in  any 
coordinate  system,  and  we  may  write 


J  ~  Bj  j  .■J'J' 
1  2'  1  2  12  1  2 


(21) 


Having  determined  the  symmetry  of  B  upon  exchange  of  initial 
and  final  states  we  now  impose  a  futher  requirement 


j;j2'-jij2 


=  BJ  J  + 
\  12 


r;j0 


/2 


(22) 


which  satisfies  (21)  and  has  the  additional  property  that  the  B  part  of 
the  T-matrix  factors  into  two  parts,  one  involving  the  initial  state  and 
the  other  involving  the  final  state.  Equation  (22)  is  convenient,  but 
not  absolutely  necessary  to  the  development  of  the  theory.  This 
somewhat  arbitrary  decomposition  will  be  shown  below,  in  any  event,  to 
affect  only  a  restricted  class  of  contributions  to  the  interruption 
function  S(b) . 


In  order  to  obtain  a  condition  on  the  A's,  we  investigate  the 
rotational  properties  of  the  T-matrix.  We  assume  that  T  may  be  expanded 
in  terms  of  unit  irreducible  tensors  in  the  spaces  of  particles  1  and  2: 


T  =  T  Tklk2  ^1(1)  ^2(2)  .  (23) 

*tqi  ^  92 

k2q2 

(These  unit  tensors  are  defined  such  that  the  reduced  matrix  elements 
between  any  states  are  unity.)  We  may  then  calculate  the  matrix 
elements  of  T: 


*12 


(24) 


<Jl("l)l‘l("l  '  MiXJ2(M2)k2(^  -  "2)|Ji  > 


And  we  may  perform  a  sum  over  the  M's  to  obtain 


<VlJ2M2|T'JlV2“2> 


12 


(25) 


‘  (2J1  *  ')(2J2  *  ’)  <V2I|T°0|IV2>  • 


This  states  the  well-known  result  that  the  average  of  the  diagonal 
matrix  elements  of  an  operator  within  the  subspace  (JjJ2)  is  e<Iual  to 
the  part  of  the  operator  which  is  invariant  under  rotations?  i.e.,  a 
scalar  in  the  separate  spaces  of  molecules  1  and  2.  We  impose  the 
constraint  that 


"i"2 


'W2VJlV2"2 


■  (2J1  *  ')(M2  +  ') 


nO 


(26) 
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where  A*n^  is  the  nth  order  term  in  the  power  series  expansion  of  A,  and 
that  the  zeroth  order  term  A^^  is  invariant  under  rotations  (is  in  fact 
equal  to  the  unit  matrix).  All  the  other  terms  in  the  expansion  of  A 
contain  no  part  which  is  invariant  under  rotations.  Referring  again  to 
equation  (12),  we  see  that  (25)  and  (26)  imply 

•  ,27) 

Therefore,  Bjjj2  is  an  exponential  expansion  of  the  rotationally 
invariant  portion  of  the  T-tnatrix  element. 

2.2  Evaluation  of  Various  Terms  in  Expansion 

We  are  now  in  a  position  to  evaluate  the  various  A*n*  and 
B^n*.  By  expanding  the  exponential  in  (12)  to  the  nth  order  and 
matching  the  kth  order  terms  obtained  thereby  with  the  kth  order  term  in 
equation  (11),  we  find 


[\k  exp(Bmk}]n  "  f  dt1  ”•  f  n"1dtn(k)n 

J-a,  J-a> 

(28) 

x  l  <m|v(tl)|i>  •••<p|v(tn)|k>  , 

**  #P/  •  •  • 

where  (ln  means  the  term  in  the  expansion  of  the  quantity  within  the 
brackets.  This  equation,  together  with  equations  (22)  and  (26),  is 
sufficient  to  determine  the  A^n^  and  B^n^  to  arbitrary  order.  We  shall 
evaluate  the  quantity  A  to  first  order  and  B  to  second  order,  and  then 
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state  a  general  prescription  for  arbitrary  order  (the  reason  for 
evaluating  A  to  one  less  order  than  B  will  become  clear  shortly). 

(0) 

We  set  B  =0.  Then  we  obtain 


.(O) 


6j .j'/jj'Vm'V 

11  2  2  11  2  2 


(29) 


as  stated  earlier.  The  first-order  term  gives 


^;J2M2’J1M1J2M2 


J'®J  J'®M  M ' ®  J ^ J 

1  1  2  2  1  1  2  2  1  1 


(30) 


Setting  { J ’ J 2M1 M2 ^  =  anc*  8unmin9  over  and  making  use 

of  equation  (26),  yields 


.M) 


J1  J2  (2J1  +  1 )  ( 2 J 2 


*  ’>  «  L 


dt. 


(31) 


„1„  <JlV2M2|,f(tl)|JlV2"2> 


12 


which,  when  substituted  back  into  (30),  gives  the  following  expression 
for  the  first-order  term  in  A. 
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and  we  may  substitute  this  expression  back  into  equation  (33)  to  obtain 
(2) 

an  expression  for  A  . 

It  is  now  clear  how  to  obtain  the  kfc^  approximation  for  A  and 
B.  We  simply  write  down  equation  (28)  for  the  kth  approximation, 

set  (J'J'M.'M')  =  (J1J2M1M2)  and  sum  over  M  and  M  .  This  gives  a 

1212  ( k)  1  L 

relation  analagous  to  (34)  for  the  B  in  terms  of  known  quantities  and 

the  B's  already  calculated.  Substitution  of  the  resulting  expression 

( k) 

into  (27)  gives  the  A  .  One  may  proceed  in  this  manner  to  any  desired 

order  in  the  expansion.  We  shall  stop  here,  however,  and  calculate  S(b) 
(k)  (k) 

in  terms  of  the  A  and  B  . 

2.3  Evaluation  of  S(b) 


We  now  substitute  our  expansion  of  the  T-matrix  into  equation 
(6)  to  obtain  an  expression  for  S(b).  To  do  this,  we  need  the  matrix 
elements  of  T  ^b),  which  are  obtained  from  those  of  T(b)  by  unitarity. 


<JiMI 


,J*M' 
1  2  2 


(b) 


<JlMlJ2M2|T(b)|JlMlJ2M2>  ‘ 


(35) 


Substituting  the  expansion  (12)  for  the  T-matrix,  we  first  consider  all 
the  terms  which  contain  a  factor  A^  .  For  these  we  may  use  the 
unitarity  of  the  V-C  coefficients  and  equation  (26)  to  perform  the  M- 
sums  explicitly.  All  the  other  terms  must  be  handled  separately,  one  at 
a  time.  With  this  in  mind,  we  arrive  at  the  following  result  for  S(b): 


19 


1  % 


1 


-  1 

J$ 

all  M 


S(b)  =  1  -  -cpp^  - 

<Jf  (Mf) 1  ( M)  1 J1  (Mi)XJf  (Mf) 1  ( M)  1  Ji  (Mi)  > 


(2Ji  +  ')(2J2  +  1_) 


[jfMf 


+  a<2>*  + 

JfMfJ2M2,JfMf J2M2  JfMfJ2M2,JfMfJ2M2 


] 


(36) 


Ld)  +  a(2)  +  1 

[VlJ^;W2M2  JiMiJ2M2,JiMiJ2M2  “*J 


x  exp[i  { 


Bt  .  +  B*  ,,  +  B- 

JfJ2  JfJ2  J 


1^2  +  V)}  ] 


Vhe  main  contribution  to  equation  (36)  arises  from  the  first  two  terms, 
the  last  term  being  a  small  correction.  In  fact,  it  is  only  the  last 
term  which  is  sensitive  to  the  decomposition  (22).  Note  also  that 
equation  (36)  is  valid  to  all  orders  of  perturbation  theory,  and  that 
contributions  to  S(b)  in  the  nth  order  come  from  the  A^k*  with  k  < 
(n  -  1),  as  noted  earlier. 


3.  EVALUATION  OF  S(b)  FOR  HARMONIC  EXPANSION  OF  POTENTIAL 

In  order  to  make  the  maxinnsn  use  of  the  formalism  we  have  developed, 
we  shall  use  the  group  theoretical  properties  of  the  interaction 
Hamiltonian.  This  Hamiltonian  may  be  written  in  terms  of  spherical 
tensors  as  follows: 


t 

< 
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V(t) 


kjk2k 

mu2 


kjk2k 


[R( t) ] 


C.  (1) 

ki^i 


C.  (2) 

k2w2 


(37) 


X  +  <kl(Wl)k2(P2)  |k(Pl  +  VZ)>  ' 

where  the  are  defined  by 


Cky(r5  =  [2k  4+  l]2  Yky(r*  '  <38) 

and  where  Ykp  is  the  ordinary  spherical  harmonic.  Equation  (37)  is  the 
most  general  form  for  the  interaction  potential  allowed  by  the  rota¬ 
tional  invariance  of  the  system;  that  is,  if  we  rotate  the  coordinates 
of  molecule  1  and  molecule  2  and  also  the  vector  R(t)  between  the  two  by 
the  same  rotation,  the  interaction  potential  must  remain  unchanged.  In 
addition,  invariance  under  coordinate  inversion  requires  that  +  k2  + 
k  be  even. 

As  an  example  of  equation  37,  we  have  the  usual  electrostatic 
multipolar  expansion,  wherein 


vklk2k[R(t)] 


(  2k  1  +  2k2  )l 
(  2k  1)  !  (2k2) 


1/2 

(”1 )k2  qk  (1)  qv  (2) 
1  K2 


(39) 


x  R(t) 


-  k1-k2-1 


k  ,kj  +  k2 


where  qk  is  a  multipolar  moment  (i.e.,  q^  is  the  dipole  moment,  q.^  is 
the  quadrupole  moment,  etc.). 


* 
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The  required  matrix  elements  of  the  interaction  potential  are  given 


by 


J2M2  |V(t)  I  J2M2> 


=  «P  [^M* 


UJ  J 

12 


)*] 


k\ K2k 
M1U2 


vkik2ktR(t)i 


x 


(2ji +  + ') 

(2jl ♦  ^"2  + f)' 


1/2 


<Jj (Olkj (0)  |JJ(0)> 


X  <j2(0)k2(0)|j-(0)  >  |j;(Mj)> 

(40) 


X<J2(M2)k2(y2)|J2(M2)><kl(yl)k2^2)|k(^  +  )  > 


x  C, 


k,^ 


(R) 

M2 


We  now  substitute  this  expression  for  the  matrix  elements  into 

equations  (31),  (32),  and  (33)  to  obtain  the  and  B^n^.  First,  we 

consider  the  8,  ,  .  The  first-order  term  (32)  is  obtained  simply  by 
J2 

summing  equation  (4)  over  the  M's: 


B 


(2) 
J1  J2 


i-fi 


at  v000  [r(  t) ) 


-iA 


(1) 
J1  J2 


(41) 
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For  the  second-order  term  we  substitute  the  matrix  elements  into  (34)  to 
obtain 


a(  2 )  _ 
BJ,J  - 


1J2  (iR) 1  (2J1  +  1)(2J2  +  1) 


+  n  L 


dt . 


JlJ2MlM2ylM2 


II1  dt2  eXp[i(WJ1J2  ~  -  fc2)]  ^^k^l) 


(42) 


x  c 


x  vk;k-k -(*2)  <J1(o)k«(o)|j»(o)><j'-(o)k1(o)|j1(o)> 

X  <J2(0)k2(0)  |J2(°)><  J2  (  0 )  k2  ( 0  )  |J2(0)> 

< J  i  (M  i)  k  l(‘p  0 1 J  i  (Mr|li)>  < J  I’  (Mi  "M  0 k  i  (u  0 1 J  i  (M  i)  > 

<  J2  ("2)k2  (-“2)  |  j;  (M2-“2)>  <■ (“2)  |J2  ("2)  > 

<kl(Ul)  R2(U2)  "'(‘,l+U2)><kl'  ('“I)  Ri(-U2),1'‘(-“l-,‘2)> 

[s(‘i)]  PM]  -  • 


k'U1+W2 


Now  all  M-summations  but  one  may  be  performed  using  the  unitarity  of  the 
V-C  coefficients  to  obtain 


b(2) 

J1J2 


/J  (0)k  (0  )  |  J"(  0  )\2  /  J  (0  )k  (0)|J"(0K2 
\  1  1  V  /  s  Z  Z  Z  / 


(2k,  +  ')(2\  *  ’ 


*£*./>.  “"{‘(V*  '  “jrcXt‘  '  t2)] 

*  c;»p(‘l)]  Ckup0=2)]  -  Bv'/2  • 


(43) 


In  order  to  perform  the  indicated  time  integrations,  we 
the  Fourier  transform  of  the  interaction  potential: 


introduce 


ki  ko 
v,  1 
kp 


(to) 


itot 

e 


V.  .  ,(t)  C. 
kjk2k  kp 


(44) 


Introducing  this  into  equation  (43)  allows  the  time  integrations  to  be 
performed.  Therefore,  equation  (43)  reduces  to 


B'“'  =  -S'  -  iA  +  —A  , 

J i  J 2  J  i J 2  JjJ2  2  J 


(45) 


where 


„  ( 2 ) outer 
3j1j2 


1  y 

2  j’i’j^p 

k  ik2k 


(0 ) kj ( 0 ) | J"(  0  )^2 
_____  _____ 


(46) 


x  <J2 (o)k2 (o)  I J2 < 0 )>2  vB«(j-j2-  -  JtJ2) 


k  l  k2  / t»»  t »» 


and 


i(2) 

J1J2 


—  I 

2lT  J\'J§P 
k  ]^k 


(JltOlkjtO)  Jj'(0)>2<J2(0)k2(0)  J2"(0)>2 
(2k  i  +  1)  (2k 2  +  1) 


X  —  P  /  dto  /u)T  T  -  u)T„Ttl  +  to ^  1  Ivjy,1,*2  (u))  I2 

^\2  J-cc 


("jj  J2  _  +  “)  V 


(47) 
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and  where  is  given  by  equation  (41): 

J1 J2 

=  J-v°°(0)  .  (48) 

J1J2  "  00 

All  the  B's  have  now  been  evaluated  to  second  order.  What  remains 
is  the  last  term  in  equation  (36),  the  so-called  correction  term.  To 
second  order  in  the  interaction  potential,  we  have 


( 2  )corr  ,  <Jf(Mf)1(M)!Ji(Mi)XJf(Mf)1<M,|Ji(Ml)> 


l 

J2’ 

all  M 


(2Ji  +  ’)  (2Jf  +  ') 


*  AWJ2M2;JfMfJ2M2  ^i«iJ2M2;JiMiJ2M2 


(49) 


expp  {BV2  +  B^fJ2  +  Bjij2  +  Bjiji } 


We  may  now  insert  equation  (31)  into  this  to  obtain 

( 2 )  corr  ,  <Jf (Mf) 1 1 Ji (Ml)> 

*2  J>  (2Ji  *  ’)(2J2  +  ’) 

2 

all  M 

X  1(M) |Ji(Mi)>  exP 


(50) 


2  i*V2  +  V*  +  V2  +  V* 


<JfMfJ2M2,9r(tl)  JfMfJ2M2>’ 


r 

‘J_.  dt2<JiMlJ2M2^(t2)|JiMiJ2M2> 


+  ^2^2  eXP[BV2  +  BJiJ2] 
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The  matrix  elements  (40)  may  be  substituted  into  the  above  and  the  M- 
summations  performed  by  standard  methods.  The  result  is 


,(2)corr  t  _(2)middle 


Ji"f '“2 


'iuf  '“2  2 


x  exp 


[i 


B*  +  B*  ,  +  B 
JfJ2  JfJ2  J 


A  *  Vi  I] 


(51) 


+  A<U  exp|~B*  T  +  B 


JiJ2  JfJ2 


■Jf  J2 


where 


(2)middle  _  _  J_  y  ^(2Ji  V)(2Jf  *  "O  w  /-jj  j  k  .j  j  \ 

Vf»J2J2  "  *2  klky  (2ki  +  ')  (ik2  +  ')  (  f  1  1' “ f) 


12 


(52) 


x  <^(0)^(0)  |Ji(0)><Jf(0)k1(0)  |Jf(0)><J2(0)k2(0)  |j£(0)>2 


and  we  have 


S(b) 


exp 


[V2 


+  B 


JiJ2J 


.  ( 2)  corr 

JiJf?J2 


(53) 


where 


B 

J1J2 


. ( 1 )  _  ( 2 ) outer 

J 1 J  2  J) J2 


iA 


(2) 

J1J2 


A(1) 

J1J2 


(54) 
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.  .  _(2)corr  .  .  „.(1)  „(2)outer 
and  where  S  is  given  by  equation  (54)  and  A  ,  S  ,  and 

(2 ) 

A  are  given  by  equations  (51),  (49),  and  (50),  respectively.  Note 
that  only  g'2^corr  sensitive  to  the  decomposition  of  equation  21;  the 
other  terms  in  equation  (53)  do  not  depend  on  it  at  all. 

It  is  appropriate  at  this  time  to  remark  on  the  similarity  between 
our  theory  and  Anderson's  theory.  Note  that  each  of  the  quantities 
appearing  in  the  above  equations  also  appears  in  Anderson’s  theory. 
They  appear,  however,  in  a  very  different  form. 


4.  THE  LINKED  CLUSTER  THEOREM 

The  following  section  is  devoted  to  the  simplification  of  the 
preceding  results  and  the  justification  of  the  exponential  expansion  by 
use  of  a  linked  cluster  theorem.  We  show  first  that  matrix  elements  of 
the  interaction  between  different  vibrational  states  may  be  ignored  as 
long  as  the  potential  contains  the  second-  and  higher-order  forces  from 
the  outset.  If  this  is  the  case,  then  the  interaction  Hamiltonian  acts 
only  on  the  rotational  variables  of  the  two  molecules.  We  then  show 
that  we  may  separate  certain  terms  from  all  orders  of  perturbation 
theory  in  the  T-opera_or,  and  that  the  terms  which  are  so  separated  may 
be  summed  to  obtain  the  exponential  expansion  obtained  earlier  by  a 
different  method. 

We  first  note  that  the  above  results  simplify  somewhat  by  observing 
that  the  Fourier  transforms  in  equation  (44)  are  negligible  when 


u  >>  2irv/b 


(55) 


<  % 
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where  v  is  the  relative  velocity  of  the  collision  and  b  is  the  collision 
impact  parameter.  Taking  v  =  3  x  104  cm/s  and  b  =  5  A,  we  obtain  a 
typical  value  of  20  cm”'  for  the  cutoff  frequency.  This  is  comparable 
to  the  rotational  splittings  of  a  molecule,  but  it  is  very  small 
compared  to  the  splitting  of  the  vibrational  levels.  Therefore, 
contributions  to  the  sums  in  equations  (52)  and  (46)  from  vibrational 
levels  other  than  the  ones  corresponding  to  ,  J  may  be  neglected. 
For  the  same  reason,  electronic  levels  other  than  the  ground  level  give 
no  contribution  to  these  sums. 

(2) 

For  the  A  term  in  equation  (47),  on  the  other  hand,  the  situation 
is  different.  For  large  values  of  the  energy  denominator,  the  integral 
in  (47)  gives  a  contribution  which  is  proportional  to  the  inverse  of  the 
energy  difference.  In  fact,  the  contribution  to  this  term  from  the  far- 
lying  vibrational  and  electronic  states  is  exactly  what  one  would  expect 
if  the  energy  of  interaction  between  the  molecules  were  calculated 
according  to  Rayleigh-Schroedinger  perturbation  theory  and  added  to  V(t) 
to  obtain  an  effective  potential.  This  summing  over  the  far-lying 
states  may  be  performed  to  arbitrary  order  in  perturbation  theory  and 
yields  the  well-known  Van  der  Waals  forces  (induction,  dispersion, 
etc.).  We  shall  assume  that  these  contributions  have  been  included  in 
the  potential  from  the  outset. 

If  we  require  that  the  higher-order  forces  are  included  in  the 
potential,  then  we  may  therefore  neglect  all  matrix  elements  of  V 
between  states  of  different  vibrational  and  electronic  quantum  numbers; 
we  therefore  consider  the  potential  to  be  a  function  only  of  the 
rotational  variables  of  the  two  molecules,  and  we  include  the  dependence 
of  matrix  elements  on  vibrational  quantum  numbers  parametrically.  Thus, 
in  a  given  vibrational  level ,  the  potential  contains  as  operators 
irreducible  tensors  which  operate  on  the  rotational  variables  of 


molecules  1  and  2.  This  reduction  of  the  vector  space  has  been  assumed 
implicitly  in  the  preceding  section;  in  particular,  it  means  that  all 
sums  over  intermediate  states  are  actually  over  states  that  belong  to 
the  same  electronic-vibrational  manifold  as  the  state  under 
consideration • 

This  reduction  of  the  vector  space  has  another  important 
consequence.  To  see  this,  we  consider  the  expansion  of  V(t)  into 
irreducible  tensors  as  given  by  (37).  We  consider  the  product  of 
several  V's: 


p(tl»t2 . fcn)  =  V(tl)V(t2)*“V(tn)  *  (56) 

Each  of  the  factors  in  (56)  is  expandable  in  the  form  of  equation  37;  we 
consider  the  products  obtained  by  taking  one  term  from  each  factor.  If 
the  tensorial  rank  of  the  first  factor  is  kj^k^  ,  the 
second  ki2*k|2'  ,  etc.,  where  kp'  is  the  rank  of  the  jth  term  in  the 
space  of  molecule  1  and  similarly  for  kpP  then  the  above  product  is  a 
sum  of  terms,  each  of  which  is  represented  by  the  Kronecker  product, 


->kl 


(1) 


x  D 


,kl 


(2) 


.x  D 


U>)  f  v(1>  1 

1  |pk2  x  D 


ik2 


(2) 


0*2 


n) 

i 


(57) 


=  l 

kLk2 


“*1*2 


Dk2 


where  the  are  irreducible  representations  of  the  rotation  group,  and 
the  coefficients  a  are  determined  by  the  properties  of  the  rotation 
group.  The  product  of  potentials  in  equation  (56)  corresponds  to  the 
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We  shall  be 


sum  of  (57)  over  all  allowable  values  of  kj^  and  k2^  • 
interested  in  separating  out  from  this  product  all  the  terms  which  have 
k)  =  k2  =  0;  i.e.,  the  rotationally  invariant  part.  Therefore,  we 
write 

. ' "(SO]"”  * 

(58) 

♦[v(t,)v(t2)...v(tn)]'  . 

where  the  first  term  corresponds  to  that  portion  of  equation  57 
with  kj  =  k2  =  0,  and  the  remainder  of  equation  (58)  corresponds  to  (57) 
with  values  different  from  kj  =  k2  =  0. 

The  first  term  of  equation  (58)  has  the  important  property  that  its 
matrix  elements  between  any  two  states  are  diagonal  in  all  quantum 
numbers.  The  fact  that  it  is  diagonal  in  the  rotational  quantum  numbers 
follows  from  its  rotational  invariance  ^kj  =  =  O^i  it  is  diagonal  in 

vibrational  and  electronic  quantum  numbers  because  of  our  reduction  of 
the  vector  space  (i.e.,  separation  of  the  higher-order  forces). 

With  these  preliminaries,  we  shall  proceed  with  the  derivation  of 
the  linked  cluster  theorem.  We  shall  define  a  linked  cluster  of  order  n 
as  a  product  of  the  form 

[?(tl)?(t2)-",(t»)]i°>  •h.OVb . ln)  •  <59> 

which  contains  no  subproducts  that  transform  like  kj,k2  =  (0*0).  *c~ 
cording  to  the  above,  a  product  of  the  form  £vf(t  ^^(*2)  *  *  ‘^^n  )]*°  ^ ' 
which  in  general  does  contain  subproducts  which  contain  kj  =  k2  =  0 
parts,  may  be  broken  up  into  parts  which  do  not  contain  any  kj  =  k2  =  0 
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parts;  i.e.,  into  a  sum  of  products  of  linked  clusters.  For  example,  we 
have  the  product 


*  [,(tl)---,(t»)]<0’  [,(lm.l)---,(tn)]<‘”  |60) 

♦[*(*1)— ' 'f.)]'  [V-.)—' ^a)]'  • 

The  first  term  corresponds  to  a  reduction  of  the  original  product  into 
two  similar  products,  which  may  themselves  be  reduced  further.  This 
procedure  may  be  continued  until  it  is  impossible  to  reduce  the  clusters 
any  further;  that  is,  until  every  part  of  equation  (60)  which  transforms 
according  to  k^  =  =  0  is  also  a  linked  cluster.  The  second  term  in 

(60)  is  not  reducible  further  as  it  stands,  but  may  be  reducible  if 
arranged  differently.  This  process  may  be  continued  until  the  entire 
product  on  the  left-hand  side  of  (60)  has  been  broken  up  into  a  series 
of  linked  clusters. 

We  now  consider  equation  (24)  for  the  symmetric  part  of  the  T- 
matrix.  We  may  substitute  this  relation  into  the  power-series  expansion 
(11)  for  the  T-matrix,  and  arrive  at  the  following  expression: 


jo  (hr 

(61) 

ft  ,  _ 

x/  dt  j  .  . .  J  n”  dtn  <J1M1J2M2|V(t1J...V(t^  |J1M1J2M2> 
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(61) 


GO  f 00  C t 


x  ...V(tn)]  <°>  UjMj  J2M2  >  . 

In  order  to  simplify  what  follows,  we  introduce  the  time-ordering 
operator  T,  which  has  the  following  property. 


T  f 


(  f('l)  f(cz)‘  h  >  *2 
(ll)  l(cz)  ■  1 


(62) 


(  '(*0  '(*>)  '  *1  <  *2 

Using  this  notation,  we  may  write  equation  (61)  as 

<JlJ2l|T°°l|JlJ2>  ■  I  f"  at,...  f"«„ 

n=0  n  v-oo  J-co 


(63) 


x  <((!», J2«2l[^(t1)...«(t„)],0lla1»1J2«2>  . 

We  now  consider  the  effect  of  decomposing  the  right-hand  side  of  this 
equation  into  linked  clusters  according  to  (60).  We  consider  the 
decomposition  of  the  nfcl1  term  into  k^  linked  clusters  of  order  1,  k2 
linked  clusters  of  order  2,  etc.  The  number  of  ways  to  do  this  is  given 
by7 


7J.  Hubbard,  Proc .  Roy.  Soc.,  A240  (1957),  539. 
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1 


k  k 
1*2 


— v - irQj - 

1  2 

(1!)  (21)  ...k  ly... 


(64) 


and  the  total  nth  order  contribution  to  (65)  is  given  by  the  sum  of  all 
such  decompositions  over  all  values  of  k^  ,  k?  ,  etc.,  subject  to  the 
restriction 


k^  +  2k2  +  3kj  +  ...=;n  . 


(65) 


When  we  perform  the  sun  over  n  we  obtain 


00  o 

klV 


dtl  [’(*■)]  "”|Jl"lJ2"2> 


(66) 


We  may  now  remove  the  time-ordering  operator,  change  back  the 
integration  limits,  and  perform  the  sum  over  the  k’s  to  obtain  the  final 
expression 


<j1j2I|t0<,||j1j2>  -  exp  [&]£  at, 


(67) 


<W2M2l[v(tl)]<°>lJ1V2>.2>e  [»]2ll  *1  £'  dt2 


<J1“lJ2M2l  [?(tl)?C2)]i'”lJ."lJ2"2  > 


+  see 
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TTius ,  we  have  rederived  the  exponential  expansion  in  a  different 
way.  If  we  compare  th_s  result  to  equation  (27),  we  see  that  the  nt^1 
order  term  in  B  is  given  by 


□  ( n) 

v2 


dt 


dtn 


(68) 


X 


1  1  2  2 


By  introducing  our  earlier  notation  into  (68)  we  see  that  is  given 

by 


,(  n) 

Va 


- (M ■ 


•  2  J . 


2J , 


E 


all  M 


(69) 


-<j;M;j2";i’(tn)|jiH.j2M2>]L  • 

where  again  the  subscript  L  indicates  that  only  linked  clusters  are  to 
be  considered;  that  is,  no  subproduct  of  the  matrix  elements  above 
contains  a  rotationally  invariant  part.  This  is  a  generalization  of  the 
Goldstone-Brueckner  linked-cluster  theorem7'®  to  degenerate  systems. 


7J.  Hubbard,  Proc.  Roy.  Soc . ,  A240  (19 57),  5 39. 
®J.  Goldstone,  Proc.  Roy.  Soc.,  A239  (1957),  267. 
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The  results  given  above  may  be  exhibited  explicitly  for  the  second- 
(2) 

order  B,  from  the  formulas  (45),  (46),  (47),  and  (48).  In 

J1 J2 

particular,  the  k  =  k  =  0  term  of  (47)  gives  no  contribution,  since 

.  1  00 
the  energy  denominator  u>T  T  u)Tn,„  vanishes  and  since  v„„(u))  is  an 

J  ^  00 
l  t  12(2) out 6  r 

even  function.  The  contribution  to  S  from  the  term  k  =  k  =  0 

2  12 

exactly  cancels  the  term  1/2A*1'  in  (45).  Therefore,  b!21  is  given  as 

12 

in  (71),  the  linked  clusters  being  precisely  the  sums  (46)  and  (47)  with 

k  =  k.  =  0  left  out.  We  shall  write  these  terms  as  g(2)outer 

(  2  )• 

and  A  respectively. 

_  ,  ....  „(2)corr 

There  is  a  further  cancellation  in  S  as  given  by  equation 

(51),  which  arises  from  the  k^  =  k?  =  0  part  of  the  sum  there. 

Therefore,  the  linked  cluster  theorem  gives  for  the  quantity  S(b)  to 
second  order. 


S(b) 


exp 


B*  - 

.  JfJ2 


+  B 


JiJ2_ 


+  S 


( 2)corr 
Ji  Jf  ;  J2 


(70) 


where 


c(2)corr  _ 
JiJf;J2  "  }, 


.  ( 2 (middle ’ 

5 Jf Jf  >j2J2 


2 


(71) 


exp  \ 


BJ  J 
JfJ2 


+  B* 


jfj; 


+  B 


JiJ2 


+  B 


JiJ' 


where  is  as  given  in  equation  (52)  with  the  k^  =  k^  =  0  part 

omitted,  and  where 


V2 


•iA(1)  - 

J1J2 


,( 2 )outer’ 

V* 


iA 


(2)* 
J  ( J  , 


(72) 
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The  linked  cluster  theorem  has  provided  both  an  alternative  way  to 
derive  the  exponential  expansion  and  also  a  simplification  of  the 
theory,  having  eliminated  extraneous  terms  from  our  equations  (45), 


(46),  (47),  (48),  and  (49).  Of  course,  we  could  have  noticed  this 

cancellation  without  going  through  the  laborious  procedure  of  deriving 
the  linked  cluster  theorem.  However,  our  derivation  has  two  advantages; 
first,  it  gives  a  much  simpler  formula  (69)  for  the  higher-order  B^n^'s 
and  second  (and  more  important),  it  provides  us  with  some  insight  into 
the  nature  of  the  exponential  expansion  and  its  relation  to  more 
familiar  expansions  of  this  type  used  in  other  branches  of  physics. 


5.  COMPARISON  TO  EARLIER  IMPACT  THEORIES 

It  is  easy  to  show  that  the  earlier  impact  theories  follow  directly 
from  the  one  presented  here,  given  certain  approximations.  The  simplest 
of  these  is  the  old  phase-shift  theory,9  which  follows  directly  from 
(70).  If  we  consider  an  intermolecular  potential  which  contains  no 
angular  dependence  whatsoever — that  is,  if  the  potential  contains  only 
k  =  k,;  =  0  in  the  decomposition  into  irreducible  tensor  operator s— -then 
our  reduction  of  the  vector  space  and  the  linked  cluster  theorem  implies 
that  all  the  B^n *  vanish  except  the  first-order  contribution,  and  all 
the  A^n*  vanish  except  A* 0 ' .  In  this  case,  we  have 

Bphase  shift  „  _iA(1>  (73) 

V 2  Jr? 


and  furthermore. 


S(b) Phase  shift  =  ,  _  e-in 


9y.  F .  ;*>■  isskopf ,  Phi/s.  Zcits,  34_  (l'>33),  1. 


(74) 


/ 
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This  is  the  familiar  phase  shift  result.9 


If  we  carry  this  type  of  analysis  further,  we  may  consider  cases  in 
which  the  potential  contains  terms  with  some  angular  dependence.  In 
this  case,  we  also  consider  the  phase  shift  separately,  but  this  time  we 
obtain  higher-order  contributions,  which  consist  of  all  the  imaginary 

parts  of  the  various  Bj  .  : 

12 


n  =  Im  B,  ,  -  Im  B, 

JfJ,  J.J2 


=  A<U  -  AiU 


J.  J 
12 


JfJ2 


A(2)  -  A(2) 

aj.j  aJ,J 

12  f  2 


(76) 


At  the  same  time  we  consider  the  other  effects  of  the  nonisotropic  part 
of  the  potential.  One  such  effect  is  the  reorientation  of  the  radiating 
molecule;  the  second  is  the  possibility  of  an  inelastic  collision.9  If 

we  single  out  all  the  terms  in  B,  T  which  have  intermediate  states  of 

J1  / 

the  radiator  which  differ  from  J  ,  then  we  may  factor  these  out  from 
S(b)  in  the  following  manner: 


Pel  =  exP 


V  Jf  j2 


)] 


(77) 


In  order  to  obtain  a  simple  expression  for  the  rest  of  *-he  terms,  we 

shall  make  the  approximation  that  in  equation  (71)  the  terms  BT  T.  and 

'  f' 

_  .  ,  _ ( 2 ) cor r 

B  T  ,  ,  respectively.  This  is  a  good  approximation  since  S  is 


9/?.  C, •  Gordon,  J.  Chem.  Phi/s .  ,  -» A  (1966),  j.  3. 
F.  Weisskopf,  Phg.s.  7.rits,  34^  (1933),  1. 
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small  to  begin  with,  and  since  the  nonresonant  terms  in  which  this 

approximation  is  made  are  small  compared  with  the  resonant  terms.  Also, 

the  B  T  .  factors  are  not  expected  to  be  strong  functions  of  J  .  With 
12  2 
this  in  mind,  the  part  of  S(b)  which  has  not  yet  been  factored  out  may 

be  written  as: 


cos2 ( a/2 ) 


exp 


[Refs' 


L  V 


JfJ2 


+  B 


'JiJ2 


)] 


(78) 


1  *  l 


c( 2 )middle ' 

Vf'Vi 


In  equation  (77),  Bl  7  represents  the  part  of  B-.  _  which  contains 

12  12 

intermediate  states  of  the  radiating  molecule  which  differ  from  J^,  and 

in  (78)  B”  ,  represents  that  part  of  B _  .  which  contains  no  such 
12  12 

states.  Clearly,  equation  (77)  represents  effects  of  inelastic 
collisions,  whereas  equation  (78)  represents  effects  of  reorientation. 

Using  equations  (76),  (77),  and  (78),  we  may  write  equation  (70)  as 

S (b)  =  1  -  P  ,  e"10  cos2 ( a/2 )  ,  (79) 

el 


which  is  of  the  same  form  as  Gordon's  semiclassical  expression  for 
S(b).  With  this  in  mind,  we  shall  interpret  as  the  probability  of 

an  elastic  collision  and  a  as  the  angle  through  which  the  angular 
momentum  is  reoriented  by  the  collision.  In  fact,  if  one  were  to  make  a 
power  series  of  equations  (76)  through  (78),  one  would  recover  Gordon's 
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equations  (42  a)  through  (42  c)  ,  in  which  Gordon's  theory  is  compared 
with  Anderson's  theory.  The  theory  presented  here  is  therefore,  in  some 
senses,  a  quantum-mechanical  version  of  Gordon's  semiclassical  theory. 
It  gives  the  three  terms  in  (79)  in  their  correct  relationship  to  one 
another;  that  is,  it  gives  factors  in  a  product  rather  than  as  three 
separate  additive  terms,  as  Anderson's  theory  does.^ 


On  the  other  hand,  one  may  obtain  Anderson's  theory  simply  by  making 
a  power-series  expansion  of  (70)  in  the  interaction  potential.  Such  an 
expansion  gives 


S(b) 


+  iA 


(  1  ) 
JfJ2 


iA 


(1  ) 
JiJ2 


A(1) 

JfJ2 


-  A(1) 

V5 


(80) 


s(2)outer'  .  stouter'  +  i4<2>  .  iA(2) 


f  2 


l  2 


1  2 


,  y  c(2)middle' 
, '  ^  J  .  J  c  :  J  J  ' 
J2  1  f'  2  2 


There  is  considerable  cancellation  in  the 
s(2)outer'  +  1  A(1)2  ,  s(2)outer  and  since 


above  result,  since 


e(2)middle'  <• 

JiJf'J2J2  "  J2J2 


A(1) 

V2 


,(1)  _  c(2)middle 


JiJ2 


“'JiJf»J2J2 


(81  ) 


Therefore,  the  effect  of  the  expansion  is  to  remove  the  primes  from  the 
(2 ) 

S  terms,  and  we  are  left  with  the  result 


S(b) 


iA(1)  - 

JiJ2 


iA 


(1  ) 
JfJ2 


+  S 


( 2 ) outer 
JiJ2 


+  S 


! 2 )outer 
JfJ2 


(82) 


,  ( 2 ) middle 

Vf'Va 


iA(2>  - 

JiJ2 


iA 


(2) 

JfJ2 


*P.  W.  Anderson,  Phys,  Rev,,  76  (1949),  647. 
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This  is  the  expression  given  by  Anderson,  modified  by  the  second-order 
phase  shift  introduced  by  Sharma  and  Caledonia. ^  The  fact  that  it  is 
derivable  from  our  theory  should  not  be  surprising,  since  both  theories 
were  obtained  from  the  same  theoretical  framework.  We  note  also  that  if 
one  further  term  were  kept  in  the  expansion,  the  phase  shifts  would  be 
multiplied  by  the  same  cutoff  used  by  Sharma  and  Caledonia  ^  for  the 
imaginary  part  of  S(b) . 

Equation  (70),  however,  is  more  than  a  simple  rearrangement  of  the 
Anderson  theory;  it  contains  a  natural  cutoff.  The  theory  of  Murphy  and 
Boggs4  also  contains  a  cutoff,  and  in  fact  the  interruption  function 
obtained  by  Murphy  and  Boggs  is  equal  to,  in  our  notation, 

s(WM»rp,y-W  .  ,  .  i  [.XP(2,..J2)  ♦ 

which  has  some  similarity  to  our  form  but  which  suffers  from  several 
defects.  First  of  all,  the  reorientation  effects  are  neglected;  i.e., 
the  2)middle  tem  is  entirely  absent  from  (83).  Second,  equation  (83) 
contains  the  effects  of  initial  and  final  states  in  an  additive  rather 
than  the  correct  multiplicative  manner.  Finally,  the  isotropic  part  of 
the  interaction  is  completely  absent  from  (83),  which  makes  the  Murphy- 
Boggs  theory  totally  inapplicable  to  rotation-vibration  spectra.  The 
last  defect  is  partially  a  result  of  the  ad  hoc  manner  in  which  M- 
averaging  is  performed  in  the  theory. 

We  have'  compared  our  theory  to  the  major  theories  of  impact 
broadening  and  shifting.  Our  theory  reduces  to  these  in  the  appropriate 
limits.  For  isotropic  interactions,  we  arrive  at  the  phase-shift 
theory.  We  have  arranged  the  terms  in  the  theory  so  that  they 


5 R .  D.  Sharma  and  G.  E .  Caledonia,  J.  Chem.  Phys.,  54_  (1970),  434. 
4J.  S.  Murphy  and  J.  E.  Boggs,  J.  Chem.  Phys.,  47_  (1967),  691 . 
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correspond  to  terms  in  Gordon's  theory.  By  expanding  out  our  results, 
we  have  arrived  at  Anderson's  theory.  We  have  shown  that  our  result  is 
similar  to  the  Murphy-Boggs  theory,  although  we  cannot  reproduce  Murphy- 
Boggs  theory,  as  we  have  reproduced  the  others.  In  the  next  section  we 
shall  compare  numerical  results  obtained  with  our  theory  with  those 
obtained  with  the  most  important  of  the  others,  the  Anderson  theory. 


6.  THEORY  AND  CALCULATIONS 


In  this  section  we  comment  on  some  aspects  of  the  theory  and  perform 
calculations.  In  particular,  we  shall  show  the  connection  between  our 
linked  cluster  theorem  and  the  Goldstone-Brueckner  theorem®  by  the 
familiar  diagrammatic  approach.  We  shall  compare  the  cutoff  obtained 
here  to  the  cutoff  introduced  by  Anderson.  We  shall  calculate  the  real 
and  imaginary  parts  of  the  broadening  cross  section  in  closed  form  for  a 
highly  idealized  case.  Finally,  we  shall  calculate  the  HC 1  self¬ 
broadening  half-widths  exactly  for  the  0  to  2  vibrational  band. 

6. 1  Diagrammatic  Representation 

It  is  possible  using  a  diagrammatic  representation  of  the 
intermolecular  interaction  to  establish  further  the  connection  between 
the  linked  cluster  theorem  obtained  here  and  the  more  familiar  theorems 
associated  with  the  many-body  problem.7'®  In  particular,  consider  the 
simple  interaction  vertex  shown  in  figure  1(a).  This  diagram  shows  two 
solid  lines  representing  molecular  states,  and  a  wavy  line  representing 
an  external  field.  Conservation  of  angular  momentum  at  each  vertex 
implies  that 


7J.  Hubbard ,  Proc .  Roy.  Soc.,  A240  (195 7),  539. 
®J.  Goldstone ,  Proc.  Roy.  Soc.,  A239  (1957),  267. 
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(84  a) 


J  +  k  =  J* 
~1  ~t  ~1 


J  +  k  =  J’  (84  b) 

~2  ~2  ~2 

This  relation  gives  the  usual  triangular  inequalities.  When 
constructing  more  complicated  diagrams,  we  must  sum  over  all  the  M- 
quantum  numbers  in  closed  loops . 

The  invariant  part  of  the  T-matrix  element  is  given  by  the  sum 

over  all  possible  closed-loop  diagrams,  as  shown  in  figure  1(b).  Note 

that  some  of  these  diagrams  are  disconnected.  The  linked  cluster 

theorem  proved  earlier  implies  the  factorization  of  such  diagrams;  the 

factorization  is  such  that  B,  ,  (see  equation  (26))  contains  no 

12 

disconnected  parts,  as  shown  in  figure  1(c). 


<  J1J2  11  T  11  >  =  1  +  (^j  +  +  ^  ^ 

(b) 


BJ1J2  = 


(c) 

Figure  1.  Diagrammatic  representation  of  intermolecular  interaction: 
(a)  basic  interaction  vertex,  in  which  solid  lines  represent  the  two- 
molecule  state.  Wavy  line  represents  external  field;  (b)  expansion 
of  isotropic  part  of  T-matrix,  which  includes  the  unlinked  graphs; 

(c)  expansion  of  B  given  by  linked  cluster  theorem,  in  which  unlinked 
graphs  are  absent. 
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6. 2  Shape  of  S(b) 

It  is  of  interest  to  investigate  the  differences  between  the 
exponential  cutoff  derived  here  and  the  straight-line  cutoff  given  by- 
Anderson.  The  interruption  functions  are  shown  for  a  typical  case  in 

figure  2.  It  can  be  seen  that  over  most  of  the  range  of  impact 
parameter,  the  exponential  cutoff  yields  a  value  for  ReS(b)  which  is 
less  than  Anderson's  ReS(b).  The  oscillation  of  S(b)  about  the  value 
unity  for  small  b  does  not  occur  in  all  cases  and  is  exaggerated  here 
for  emphasis.  In  general,  it  is  a  very  small  effect  unless  the  phase- 
shift  terms  completely  overwhelm  the  gt  pouter  terms. 


Figure  2.  Shape  of  interruption 
function:  curve  (a)  this  work, 

curve  (b)  Anderson  theory. 


1P.  W.  Anderson,  Phys.  Rev.,  76_  (1949),  647. 

^P.  W.  Anderson,  PhD  Dissertation,  Harvard  University  (1949). 
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6.3  Evaluation  of  Cross  Sections  for  Special  Cases 

Under  certain  conditions,  the  integral  of  S(b)  over  the  impact 
parameter  may  be  performed  in  closed  form.  This  will  be  true  in  general 
if 


(1)  A  given  term  in  the  potential  dominates  all  other  terms, 

and 

(2)  For  these  terms,  the  only  important  contributions  to  S(b) 
are  the  Stark  and  resonance  terms. 

Under  these  conditions,  the  only  nonzero  factors  in  S(b)  are  A^1^,  the 
first-order  phase  shift,  and  g^)outer^  the  second-order  broadening 
term.  Each  of  these  terms  will  be  proportional  to  some  power  of  the 
impact  parameter. 

The  largest  terms  in  the  potential  which  contribute  to 
A^1^  are  the  induction  and  dispersion  terms,  which  are  proportional  to 
the  inverse  sixth  power  of  the  intermolecular  distance.  We  mav  write 


V000(R)  =  -A/R6 


(85) 


where  the  constant  A  depends  on  the  vibrational  quantum  numbers.  The 
integral  (41)  for  the  first-order  phase  shift  may  be  readily  evaluated 
for  this  potential  to  give 


A(1)  =  -3TTA/(8fivb5) 


(86) 
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To  evaluate  g(2)outer  ^  consider  as  an  example  the  dipole- 
dipole  interaction  between  two  linear  molecules.  We  may  write  this 
interaction  in  terms  of  spherical  tensors  as 


/6MlM2 


'd-d 


I  C,m  <2> 


,3  mjm2  imf  lm2 


x  cl 


'2  *m  1+m2  (R>  <1(mi)1(m2)|2(mi  +  m2» 
and  therefore  we  have,  from  equation  (36), 


(87) 


VU2(R)  =  /6  MtM2/R3  .  (88) 

Now  the  summation  may  be  performed  in  (46)  using  this  potential,  and  the 
integrals  in  (44)  may  be  evaluated  to  give6 

smouter  <Jj  (0)K0)  lj;(0)>2 

(89) 


x  ^J2(0)1(0) |J^(0)>2  f  1(k>  ^  , 

where 

f^k)  =  £  k4  [x2(k)  +  4K2 (k)  +  3K2(k)J  (90) 


6C.  J.  Tsao  and  B.  Curnutte,  J.  Quant.  Spect.  Rad.  Trans.,  2_  (1962), 

41. 
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1.  m 


and  where 


k  =  (wj;j«  -  <‘>J1J2)b/v 


The  V-C  coefficients  vanish  for  Stark  splitting  =  J'^,  J?  =  J1^.  Iff 
however,  the  molecules  are  identical,  then  one  may  obtain  k  =  0  by 
setting  J2  =  J’  and  J2  =  (resonance).  We  shall  assume  that  all  other 
contributions  are  small  compared  to  the  resonance  contribution.  In  this 
case,  since  f^  (k)  is  unity  for  k  *  0,  we  obtain: 

/u  i  W  2\  2 

SjJj2UtSr  =  IV^/  ^7  <J1(0)1(0)  ,J2(0)>2<J2(0)1(0)  |J1(0)>?'  '  (92) 


where  J,,  =  ±  1.  Since  the  second-order  phase  shift  vanishes  on 

resonance,  and  since  the  g^middle  term  vanishes  for  dipole-dipole 
interactions,  we  have  determined  all  contributions  to  S(b)  and  we  may 


S ( b)  =  1  -  exp£-i(a0/b)5  -  (bQ/b) 


where  aQ  and  bQ  are  given  in  terms  of  the  coefficients  in  A  and 

<2)outer 
S  . 


Using  the  expression  (92)  in  equation  (4),  the  resulting 
integral  may  not  be  evaluated  in  closed  form.  However,  if  aQ  <<  bQ,  we 
may  expand  the  imaginary  part  of  the  exponential  to  obtain 


a  ~  2n 


jp,  J.  -  .-(b”/b)4  *  1(a/b)^(Vb)“  . 
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Evaluation  of  the  integral  yields 


o  = 


=  2ub02 


(95) 


On  the  other  hand,  we  may  use  the  Sharma-Caledonia^  cutoff  in  the 
Anderson  theory  to  obtain 


a.  =  2  it 
A 


j[  bJb  j  (b0/bV 


+  i 


i(a0/b)' 


[’  -(b»/b)'] 


(96) 


which  may  be  evaluated  to  yield 


a 


A 


=  2irbQ 


=  2trb02  1  +  i  |j(a0/b0)5 


2  1  +  1(0.1905)  (ao/bq)5 


(97) 


Thus,  there  is  a  10-percent  difference  between  the  real  parts  of  the 
cross  section  in  the  two  theories,  whereas  Anderson's  imaginary  part  is 
a  factor  of  1.6  lower  than  ours.  If  we  consider  velocity  averaging, 
this  will  increase  the  imaginary  part  of  the  cross  section  by  a  factor 
of  1.33  in  both  theories,  but  will  have  no  effect  on  the  real  part  of 
the  cross  section . 
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!  6.4  Exact  Calculation  for  HCL  Self-Broadening 

As  an  illustration  of  the  theory,  we  have  performed  exact 
calculations  of  the  self-broadening  in  the  0  to  2  vibrational  band  of 
HC1.  Interactions  considered  in  the  calculation  were  the  dipole-dipole, 
dipole-quadrupole ,  quadrupole-dipole ,  quadrupole-quadrupole ,  dipole- 

octupole,  and  octupole-dipole  interactions.  Also  considered  were  the 
isotropic  parts  of  the  induction  and  dispersion  interactions,  for  which 
the  difference  of  coefficients  A  in  -A/r®  was  taken  to  be 
12280  cm  1  A®.  Other  contributions  to  the  interatomic  potential  (which 
may  be  important  for  foreign  gas  broadening  or  for  cases  where  there  is 
I  no  permanent  dipole  moment)  were  neglected. 

Contributions  to  S(b)  from  the  various  interactions  were 
computed  as  in  Tsao  and  Curnutte^  and  Isnard  et  al.*®  Equation  (79)  was 
used  for  S(b)  .  Molecular  parameters  used  in  the  calculation  cure  those 
given  in  table  1.  The  quantities  u,  0,  0,  and  B  are,  respectively,  the 
(  dipole  moment,  quadrupole  moment,  octupole  moment,  and  rotational  energy 

(constant.  Results  of  the  calculation  are  shown  in  table  2  and  figure 

3.  In  the  table,  |m|  is  equal  to  max^,J^  according  to  the  usual 
notation.  We  also  show  the  results  of  a  calculation  using  the  same 
j  parameters  and  using  the  Anderson  cutoff.  The  experimental  numbers  are 

an  average  of  the  data  giver,  in  Toth  et  al11  and  Smith  et  al.12 


6C.  J.  Tsao  and  B.  Curnutte,  J.  Quant.  Spect.  Rad.  Trans.,  2_  (1962), 

41. 

10P.  Isnard,  C.  Boulet,  and  A.  Levy,  J.  Quant.  Spec t.  Rad.  Trans., 
(1973),  1433. 

llR.  A.  Toth,  R.  *.  Hunt,  and  E.  K.  Plyer ,  J.  Mol.  Spect.,  35_  (1970), 
110. 

12A.  Levy,  E.  Pioljet,  J .  P.  Bouanich,  and  C.  Haeusler,  J.  Quant. 
Spect.  Rad.  Trans.,  10  (1970),  203. 
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TABLE  1. 

MOLEClfLAR  PARAMETERS  IJSEn 

IN  THE  CALCULATIONS 

TABLE  2.  COMPARISON  OF  EXPERIMENTAL 

THEORETICAL  RESULTS  Row  HOI 

BROADENING 

Units  for  widths  are  lo”^  cm”1 /atm. 

AND 

SELF- 

Quan  t i ty 

Upper* 

state 

Lower* 

state 

|  Calculation* 

u  ( n ) 

1.  167 

1.  108 

|m| 

Exper iment  *  1  '  1  ? 

This 

AnHer  son 

e  (D-i) 

3.93 

3.60 

work 

t  heory 

i!  (D-A2) 

4.0 

4.0 

1 

21S 

mm 

2 

226 

B  ( cm”  ‘  ) 

9.82ft 

10.440 

3 

234 

mm.. 

•Values 

of  the  dipole  moment 

4 

226 

223 

242 

eere  taken  from  Smith.1*  The 

5 

207 

201 

217 

equil ibrium  value 

for  the 

*Cal 

cul at o< i  usinj  the 

interactions 

loser  ibt 

guadrupole 

moment  is 

&?Leeuw  and 

in  the 

text .  Molecular 

parameters  are  t/iverj 

Dynamus , 1 4  and  the  first  term  in 
its  vibrational  dependence  is 
from  Friedr.unn  and  Kinel.'L'>  The 
value  of  the  octupcle  moment  is 
taken  from  Sharna  and  Caledonia 3 
it  is  just  an  est imate ,  but  seems 
much  more  reasonable  than  that  of 
Isnard  et  a/.10  Rotational 

constants  are  from  Smith. 13 


table  1. 


Figure  3.  Self-broadened 
halfwidths  of  HCl:  Solid 
line,  this  work.  Dashed 
line,  Anderson  cutoff. 
Circled  points,  experimental 
data  of  Toth  et  al^ 1  and 
Levy  et  al . 12 
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Agreement  between  our  theory  and  the  experimental  data  is 
excellent;  moreover,  it  is  much  better  than  the  agreement  between 
Anderson's  theory  tnd  experiment.  In  earlier  calculations10  using  the 
Anderson  theory,  it  was  necessary  to  invoke  an  abnormally  large  value  of 
the  octupole  moment  (15  D-A^)  in  order  to  match  the  experimental 
results . 


7.  CONCLUSIONS 

A  theory  of  impact  broadening  and  shift  has  been  developed  which 
combines  the  best  features  of  the  Anderson,  Gordon,  and  Murphy-Bogqs 
theories  in  that  it  uses  Anderson's  expression  for  S(b)  (equation  (fell, 
which  is  exact  in  the  impact  limit  and  employs  an  exponential  expansion 
of  the  T-matrix  to  ensure  convergence  of  integrals  over  the  impact 
parameter.  The  final  form  for  S(b)  contains  only  terms  which  are  in  the 
Anderson  theory,  but  these  terms  are  arranged  in  such  a  way  that  the 
final  form  (79)  for  S(b)  is  similar  to  that  of  Gordon.  Evaluation  of 
the  broadening  and  shift  cross  sections  for  certain  simple  cases  shows 
that  whereas  the  broadening  is  not  particularly  sensitive  to  the  cutoff 
procedure,  the  shift  is.  Similar  conclusions  apply  to  the  influence  of 
velocity  averaging.  Calculation  of  broadening  parameters  exactly  for 
HCl  self-broadening  indicates  significant  improvement  over  calculations 
performed  using  Anderson’s  theory. 
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